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ABSTRACT 

For decades, cosmologists have been using galaxies to trace the large-scale distribution of 
matter. At present, the largest source of systematic uncertainty in this analysis is the challenge 
of modeling the complex relationship between galaxy redshift and the distribution of dark 
matter. If all galaxies sat in the centers of halos, there would be minimal Finger-of-God (FoG) 
effects and a simple relationship between the galaxy and matter distributions. However, many 
galaxies, even some of the luminous red galaxies (LRGs), do not lie in the centers of halos. 
Because the galaxy-galaxy lensing is also sensitive to the off-centered galaxies, we show that 
we can use the lensing measurements to determine the amplitude of this effect and to deter- 
mine the expected amplitude of FoG effects. We develop an approach for using the lensing 
data to model how the FoG suppresses the power spectrum amplitudes and show that the cur- 
rent data implies a 30% suppression at wavenumber k = 0.2 hMpc^ 1 . Our analysis implies 
that it is important to complement a spectroscopic survey with an imaging survey with suffi- 
cient depth and wide field coverage. Joint imaging and spectroscopic surveys allow a robust, 
unbiased use of the power spectrum amplitude information: it improves the marginalized error 
of growth rate f g = d In D/d In a by up to a factor of 2 over a wide range of redshifts z < 1.4. 
We also find that the dark energy equation-of-state parameter, wq, and the neutrino mass, /„, 
can be unbiasedly constrained by combining the lensing information, with an improvement of 
10-25% compared to a spectroscopic survey without lensing calibration. 
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1 INTRODUCTION 

Over the past three d ecades, astronomers have been conducting ever l arger redshift surveys in their efforts to probe the large-sca le struc 



ture of the universe dDavis & Huchralll982t Ide Lapparent et au[l986l : iKirshner et all 1 19871; I York et al. Eoool : iPeacock et al.l[2001h . In the 
coming decade, we a re embarking on even larger surveys: BOSsFi Wiggle j 2 \ telake et alj|201ll) . ViperjH FMOsFl HETDEjfl. BigBOSS^l 



l Schlegel et ai1l2009h . LAMOSTQ Subaru PFS-fl Euclid^ and WFIRSlJiS This upcoming generation of surveys are motivated by our desire 
to understand cosmic acceleration and to measure the composition of the universe by simultaneously measuring geometry and dynamics. The 
combination of cosmic microwave background (CMB) data and large redshift surveys trace the growth o f structure formation from the last- 
scattering surface (z ~ 1100) to low redshifts and determine cosmological parameters to high precision dWang et al ] | 19991 : lEisenstein et al 



1 http://cosmology.lbl.gov/BOSS/ 

2 http://wigglez.swin.edu.au/site/ 

3 http://vipers.inaf.it/ 

4 http://www.naoj.org/Observing/Instruments/FMOS/ 

5 http://hetdex.org/ 

6 http://bigboss.lbl.gov/ 

7 http://www.lamost.org/website/en 

8 http://sumire.ipmu.jp/en/ 

9 http://sci.esa.int/euclid 

10 http://wfirst.gsfc.nasa.gov/ 
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19991: iTegmark et ai1l20o4 ICole et aljliool) . Measurements of the baryon aco ustic oscillation (BAO) scale provide us w ith a robust geo- 
metrical probe of the angular diameter distance and the Hubble expansion rate (Eisenstein et al Perci val et al.ll2007h. Observations of 
redshift-space dis t ortion measure the growth rate of structure formation Jzhang et al.l2007 ; Guzzo et al.l2008l:|Wang|20 08; Guzik et alj201ol : 



White et alJ2009t IPercival & White1l2009l : ISong & Perciva1l2009l : ISong & Kavcll20ld : 



Yamamoto et al. 



20ld : lTang etal 



201 If ). Combining 



measurements of the growth of structure formation and the geometry of the universe prov ides a key clue to understanding the natur e of dark 
energy, properties of gravity on cosmological scales, or the nature of cosmic acceleration (lAlbrecht et al . 2006; Pe acock et al. 2006). 

The galaxy power spectrum in redshift space, a dire ct observable from a redshift survey, is a two-dimensional fun ction of wavelengths 
perpendicular and parallel to the line-of-sight direction jPeacock et al JI200 ll : lOkumura et a l. 2008; Guzzo et al.ll2008l) . While galaxy clus- 
tering in real space is statistically isotropic in an isotropic and homogeneous u niverse, the line- of-sight components of galaxies' peculiar 
velocities alter galaxy clustering in redshift space dKaiserlll987t). For review, s ee Hamilton ( 1998). The amplitude of the distortion depends 



both on geometry and dynamics I Alcock & Paczvnskill 19791 ; ISeo & Eisensteirj|2003l) . 



For the surveys to achieve their ambitious goals for precision cosmology, we will need a detailed understanding of the underlying 
systematics. One of the major systematic uncertainties in redshift-space power spectrum me asurements is non-linear redshift d istortion due 
to the internal motion of galaxies within halos, the so-called Finger-of-God (FoG) effect Jjackson 1972 ; Scoccimarrol 2004 ). Since it is 
sensitive to highly non-linear physics as well as difficult to model galaxy formation/assembly histories, the FoG effect is the dominant 

syste matic in redshift surveys. 

Reid et al. 1 20091) advocated using halos rather than Luminous Red Galaxies (LRGs; Eisenstein et alj2001 ) to trace large-scale structure. 



In an analysis of LRGs sampled with the Sloan Digital Sky Survey (SDSsJ^ iReid et alT i 20 id) implemented this scheme b y removing 



satelli te LRGs from the same halo with the aid of the mock catalog and the halo model prescription. From the SDSS LRG dataset, Reid et all 
(201(3) ound that about 6% of LRGs are satellite galaxies, while the remaining 94% are central galaxies of halos with masses 



10 13 M< 



©• 



Once such a halo catalog is constructed, clustering properties of halos are easier to model, because halos have only bulk motions in large-scale 
structure, and therefore have the reduced FoG effect. Despite this effort, the remaining FoG effect is a dominant systematic uncertainty. 

FoG effects are just one of the non-linear systematics. Future analysis of the redshift-space power spectrum of halos will need to 



refined perturbation theory suggest that halo clustering based approach seems a 


very promising probe of cosmology i 


Scoccimarrol 20041 


Crocce & Scoccimarro 2006; Matsubara 2008; Saito et al. 201 1 ; Taruva et al. 20K 


i Tana et al.l201 1 ; Reid & Whitehoi 1 


: Sato & Matsubaral 


201 lb. 



For a halo-based catalog, a significant source of uncertainty is the position of the galaxies in the halos. IHq et alj |2009) compared 
LRG positions with the X-ray su rface-brightness peak, reporting a sizable positional difference. For the LRG analysis, this is the dominant 
uncertainty (see Reid et ai]l201oL for the useful discussion in Appendix C). The virial theorem implies that off-centered LRGs are moving 
relative to the halo center thus producing an FoG effect. 

In this paper, we propose a novel method of using a cross-correlation of spectroscopic galaxies (e.g., LRGs) with background galaxy 
images to correct the FoG contamination to the redshift-space power spectrum. Dark matter halos hosting spectroscopic galaxies induce a 
coherent lensing distortion effect on background galaxy images, and the signals are measurable using the cross-correlation method - the so- 
called galaxy-galaxy or cluster-galaxy weak lensing. The lensing signals have been now measured at a high significance by various groups 
(Mandelbaum et al. f l2006l ; Isheldon et alj|2009l ; ^eauthaud et al. I l2010l : bkabe et aljboid) . S we include off-centered galaxies and use the 
galaxy position as a h alo center proxy of eac h halo in the lensing analysis, the lensing signals at angular scales smaller than the typical offset 
scale are diluted (see lOguri & Ta kada 2011, for a useful formulation of the off-centering effect on c luster-galaxy weak lensing). Thus the 



Leauthaud et al 



galaxy-galaxy lensing si gnals can be used to infer the amount of the off-centered galaxy contamination ( Johnston et al .120071 : 
20ld;iokabe et alJ20ld) . Furthermore, since lensing is a unique means of reconstructing the dark matter dis tribution, it may allow us to infer 
the halo center on individual halo basis if a sufficiently high signal-to-noise ratio is available (O guri et al-lEoioh . Hence, a weak-lensing 
based calibration of the FoG effect in redshift-space power spectrum measurements may be feasible if spectroscopic and imaging surveys 
observe the same region of th e sky. Fortunately, ma ny upcoming surveys will survey the same region of the sky: the BOSS and Subaru Hyper 
SuprimeCam (HSC) Survey (Mivaza ki et al . 2006), the Subaru P FS and HSC surveys (Subaru Measureme nts of Images and Redshifts: the 
SuMIRe project), Euclid and WFIRST or a combination of LSST 1 LSST Science Collaborations et al.ll2009h with spectroscopic surveys. 

In Section[2] we will first develop a model of computing th e redshift-spa ce p ower sp e ctrum of LRGs based on the halo model approach 
(see Coorav & Shetbj 2002. for a thorough review). Extending Whitejj 2001 ) and Seliak 2001), we model the distribution of off-centered 
LRGs as a source of FoG distortions. Following the method in lO guri & Tak ada (201 1), we also model the distribution of off-centered LRGs 
as a source of smoothing of the LRG-galaxy lensing signal. Assuming survey parameters of the Subaru HSC imaging survey combined with 
the BOSS and/or Subaru PFS spectroscopic surveys as well as the Euclid imaging and spectroscopic surveys, we study the impact of the 
FoG effect on parameter estimations. We also study the ability of the combined imaging and spectroscopic surveys for correcting for the 
FoG effect contamination based on the off-centering information inferred from the LRG-galaxy lensing measurements. For the parameter 
forecast, we pay particular attention to the dark energy equation-of-state parameter, wo, the neutrino mass parameter, /„, and the growth 
rate at each redshift slice. Unless explicitly stated we will throughout this paper assume a WMAP-normalized ACDM model as our fiducial 



http://www.sdss.org/ 
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Figure 1. A schematic illustration of the redshift-distortion effect on redshift-space clustering of dominant luminous red galaxies (DLRGs; see text for details). 
We assume that a catalog of DLRGs is constructed so that each halo contains one DLRG. The redshift distortion effect on the redshift-space power spectrum of 
DLRGs arises from two contributions: the bulk motion of halos that host each DLRG, and the internal motion of DLRG within a halo, the Finger-of-God (FoG) 
effect. If some of the LRGs are not in the center of their halos, then their motions produce significant FoG effects. The halo bulk motion causes a displacement 
of halo position in redshift space, while the internal motion stretches the distribution region of DLRGs within a halo along the line-of-sight direction. 



cosmological model ItComatsu et aijl 2009): U h h 2 = 0.0226, V. cd mh 2 = 0.1109, fi A = 0.734, respectively, r = 0.088, n a = 0.963, 



A(k = 0.002Mpc -1 ) = 2.43 x 10~ 9 , where fi b , flcdm and Q A are the energy density parameters of baryon, CDM and dark energy (the 
cosmological constant with too = —1 here), r is the optical depth to the last scattering surface, and n 3 and A are the tilt and amplitude of 
the primordial curvature power spectrum. 



2 FORMULATION: REDSHIFT-SPACE POWER SPECTRUM 



In this section, we give a formulation for m odeling the redshift-space power spectrum of luminous red galaxies (LRGs) based on the halo 
model approach JwhitehoOlLISeljakhoOll) . 



2.1 Dominant Luminous Red Galaxies (DLRGs) 



2 fi 



Wake et al. 200 



Weak lensing studies (Mandelbaum et al. 2006; Johnston et al. 2007) and clustering analyses (Ross et al. 2007, 2008. 
Izheng et al.ll20Qgl : iReid & Spergelll2009l : lwhite et al.ll201lh find that most LRGs reside in massive halos. While the typical massive halo 
contains only one LRG, roughly 5-10% of all LRGs are satellite galaxies in a halo containing multiple LRGs. These satellite galaxies 
contribute a large one halo term that is an additional source of shot noise and non-linearity in power spectrum estimation. Reid et al. 1 20091) 
outline a procedure of identifying these satellite LRGs through finding multiple pairs that lie in common halos (or the small spatial region) 
and then using only the brightest luminous red galaxies in each halo as a tracer. We call thes e galaxies domina nt l uminous red galaxies 
(DLRGs). These DLRGs are more linear tracers of the underlying matter field than the LRGs are. Rei d et al. I fcoid) and lPercival etafl bold) 
adopt this procedure to determine the SDSS LRG power spectrum. In this paper, we focus on these DLRGs so that each halo contains either 
zero or one DLRG. 



2.2 Halo Model Approach for DLRGs 



Since there is only one DLR G per halo, the two-halo term determines the clustering of these galaxies in the halo model picture l lCoorav & Sheth 
l2002l :l Takada & J ain 2003). If the DLRGs sat in the center of each halo, then the DLRG power spectrum would be linearly related to the halo 
power spectrum. However, since the DLRGs do not always lie in the center of the halo, the power spectrum is given as 

PdlrgW = f dM I dM 1 ^N H oo(M)p oS (k; M)-^Nnon(M')p oS (k; M')P hh (k; M, M'), (1) 

D L R.G J J 

where dn/dM is the halo mass function, Nhod{M) is the halo occupation number (note Nhod ^ 1 as described below), and Phh(k; M, M') 
is the cross-power spectrum of halos of masses M and M' . Numerical simulations show that the halo cross-power spectrum is approximately 
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a linearly biased version of the matter power spectrum dReid et al 1 l2009h : P hh (k; M, M') ~ b{M)b(M')P™ h (k), where 6(M) is the halo 



bias, and P 1 „ L (fc) is the non-linear matter power spectrum. This approximation simplifies the relationship between the DLRG and matter 
power spectrum: 

-I 2 



dM ^fb{M)Nnou{M)p oS {k;M) 



P£ L (fc). (2) 



riDLRG J dM 
The quantity Sdlrg is the mean number density of DLRGs defined as 

"DLRG = j dM ^Nhou(M). (3) 
The mean bias of halos hosting DLRGs is defined as 

b = — - — / dM b(M) 4r? N noo (M) . (4) 

W.DLRG J dM 

The mean mass of halos hosting DLRGs is similarly estimated as Mh = (1/?i-dlrg) / dM M(dn/dM)NuoD(M). 

The coefficient p s(k; M) in Eq. l|2} is the Fourier transform of the average radial profile of DLRGs within a halo with mass M: 

p oS {k;M) =4tt / r 2 drp oB (r)^j^-, (5) 
Jo kT 

where p g(r) is normalized so as to satisfy J vlr 4irr 2 dr p s(r) — 1, and r v i r is the virial radius of a halo with mass M, which can be 
defined once the virial overdensity and the background cosmology are specified. Note that, since the power spectrum is a statistical quantity, 
we just need the averaged DLRG distribution within a halo, which is therefore a one-dimensional function of radius r with respect to the halo 
center in a statistically homogeneous and isotropic universe. In the following, quantities with tilde symbol denote their Fourier-transformed 
coefficients for our notational convention. 

The term in the square bracket in Eq. (0 describes the halo exclusion effect. Because halos have finite sizes, roughly their virial radius, 
there is only one dominant galaxy in this region (see Fig. [T}. If we are implementing an algorithm that eliminates multiple galaxies in the 
finite region, then we impose an exclusion region around each galaxy. The two halo term describes the correlations between two DLRGs 
in two different halos. The DLR G power spectrum at s mall scales (large k's) is thus suppressed compared to the matter power spectrum 
multiplied with b 2 (see Fig. 1 1 in lCoorav & Shetbll2002h . 



If each DLRG resides at the center of each halo (e.g., the center of mass), p q(k) = 1 (or p fl oc 5n{r)). However some fraction of 
DLRGs in the sample are expected to have an offset from the halo center ( Skibb a et al l|). Due to the collision-less nature of dark matter, 



dark matter halos lack clear boundary with surrounding structures and do not have a spherically symmetric mass distribution. Thus the halo 
center is not a well-defined quantity. While DLRGs, the most massive galaxy in the halo, will eventually sink toward the center of the halo 
through dynamical friction, many clusters are dynamically young and have experienced recent interactions. Thus, we expect that DLRGs are 
not all in the centers of halos and that the distribution of the their positions in the halos evolve with redshift. 

How does this halo model picture need to be changed in redshift space? To model the redshift-space power spectrum, we need to properly 
take into account the redshift distortion effect due to peculiar velocities of DLRGs. If all DLRGs are located at the center in their host halos, 
DLRGs move together with their host halos having coherent, bulk velocities in large-scale structure, and the redshift-space clustering is not 
affected by the FoG effect. However, as illustrated in Fig.Q] if some DLRGs are offset from the center, they will have internal motions within 
their host halos, which causes the FoG effect. The virial theorem implies that the amplitude of the displacement of the DLRG from the center 
of its halo is directly related to the DLRG velocity dispersion within the halo. 

In the halo model picture, the FoG effect can be incorporated by stretching the average radial profile of DLRGs along the line-of-sight 
direction by the amount of the internal motion, as illustrated in the right panel of Fig. [TJ This stretch enhances the halo exclusion effect, 
which suppresses the power spectrum amplitudes. Thus the redshift-space distribution of DLRGs within a halo becomes two-dimensional, 
given as a function of two radii, r± and r«, perpendicular and parallel to the line-of-sight direction with respect to the halo center. Also note 
that the internal velocity distribution of DLRGs within a halo is considered to be statistically isotropic and therefore it depends on the radius 
r from the halo center, halo mass M and redshift z (see Section l2.3.3l for details). The averaged redshift-space distribution of DLRGs within 
a halo, denoted as p s ,off(f_L, fii), can be given as a smearing of the real-space distribution with the displacement function: 

/oo 
dr\ ii(ru - rf| ; r', M) p oS Ur 2 ± + rf, 2 J , (6) 
- oc 

where iJ(Arii; r, M) is the displacement function of DLRGs due to the velocity distribution inside a halo and satisfies the normalization 
condition: Jd(Ary ) R(Ar\\ ) = 1. Assuming that the internal motion of DLRGs is much smaller than the speed of light, the displacement 
of the radial position of a given DLRG is directly related to the line-of-sight component of the internal velocity dii as 

= J(ry (7) 

where H(z) is the Hubble expansion rate at the redshift of DLRG. 
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Hence, assuming a distant observer approximation, the redshift-space power spectrum of DLRGs can be given in terms of the Fourier- 
transform of p at0 f( (ri, rn ; M) as 

1 2 



P s ,dlrg(&, m) = 



"DLRG 



flu 

dM _6(M)iV H oD(M)p 3 , off (fc, fi; M) 



(8) 



where fi is the cosine angle between the line-of-sight direction and the wavevector k, i.e. /i = fc||/fc, and P^(k±,k\\) is the non-linear 
redshift-space power spect rum. In this p aper, we simply assume that the redshift distortion effect due to the coherent bulk motion of halos is 
described by linear theory ( Kaiseill987 ): 



(9) 



where /3 = f g /b, f g is the linear growth rate, f g = dlnD/dlna, and b is the effective bias of halos hosting the DLRGs (Eq. |4)). As given 
by the term in the square bracket in Eq. (|8j, the FoG effect due to off-centered DLRGs causes scale-dependent, angular anisotropies in the 
redshift power spectrum amplitudes. 

In the limit that all DLRGs are at the true center of each halo, the redshift-space power spectrum (Eq. JH]) is reduced to the halo power 
spectrum in redshift space: 

I 



Poff(r) 



Airr 2 



S D (r) 



s.DLRG 



(k,n) = b P S Jn (k,n) 



(10) 



where P s ,haio(fc, is the halo power spectrum in redshift space. More rigorously speaking, halo clustering is affected by non-linearities in 
gravitational clustering, redshift distortion and biasing at scales even in the weakly non-linear regime we ar e interested in IScoccimarrol2004l 
Taruv a et al.l2009ll20ldlSaito et alj|201 ll ; lTang et al.l201 ltlReid & Whitel201 ll ; ISato & Matsubar3l201 ll) . For these we can use an accurate 



mode l of the redshift-space sp ectrum of halos by using refined perturbation theory and/or N-body and mock simulations (ITaruva et al .120091 
20 id : Sato & Matsubaralboi lh . Hence we can extend the formulation above in order to include these non-linear effects, simply by using 



a model of non-linear, redshift-space halo power spectrum for P s ,hh(k, fx; M, M'), instead of b(M)b(M')P^^(k, fi) in Eq. (8). However, 
this is beyond the scope of this paper, and we here focus on the FoG effect by assuming the Kaiser formula (|9) for the sake of clarity of our 
discussion. 



2.3 Model ingredients 

To compute the redshift-space power spectrum (Eq. JH]), we need to specify the model ingredients: halo occupation distribution of DLRGs, 
the off-centered distribution of DLRGs and the velocity distribution inside halos. In this subsection, we will give these model ingredients 
adopted in this paper. 



2.3.1 



HOD and the halo model ingredients 

First we need to specify the halo mass function and the halo bias. We use the fitting fo rmula developed in Sheth & Tormen ( 1999h to compute 
the halo mass function and the halo bias in our fiducial cosmological model (also see iTakada & Jamll2003h . 



A useful, empirical method for describing clustering properties of galaxies is the halo occupation distribution (HOD) dScoccimarro et al 



2001; Zheng et al. 2005, also see references therein). The HOD gives the average number of galaxies residing in halos of mass M and at 
redshift z. The previous works have shown that the halo model prediction using th e HOD modeling can well reproduce the observed propertie s 
of LRG clustering over wide ranges of length scales and redshifts (0 £ z £ 0.5) l Zheng et alj2009l : lReid & Spergell2009l:lwhite et alj201 lh . 



Since w e assume that sa t ellite LRGs can be removed based on the method of iReid et al 
found in Reid & SpergellJ2009h : 



i2Qld), 



iV HO D(M) ~ N cen (M) 



1 + erf 



log 10 (Af)-log 10 (M n 

O log M 



we use the HOD for central LRGs that is 



(11) 



where erf (a;) is the error function, and we adopt M m i n = 8.05 x 1O 13 M0 and <ji og a/ = 0.7. We do not consider a possible redshift evolution 
of the HOD, because any strong redshift dependence has not been found from actual data. Note A r HOD(Af) ^ 1. Also note that we use the 
HOD model for "central" galaxies, but this does not mean that all DLRGs under consideration are central galaxies, but each halo has one 
DLRG at most. 



2.3.2 Radial profile of DLRGs 



The radial profile of DLRGs is not well known, as the true center of a halo is not easy to estimate observationally. Several studies, both 
observational and numeric al, suggest that the DLRG are more centrally concentrated than the dark matter, but do not all lie in the bottom of 
the d a rk matter potent ials (Lin & Mohr 



2004 



Koester et al 



2001 Ijohnston et al]|2007l ; Iho et aljlioolil : Iffilbert & Whitelboioj : bkabe et al 



20 id lOguri et alJlioiol ; Iskibba et alJboilh . Iho et alj feoCCT compared the LRG positions with X-ray peak positions for known X-ray 
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Figure 2. This figure shows the statistically averaged radial profile of dark matter and DLRGs in halos with mass M = 1O 14 /i -1 M0 and halo concentration 
c vir = 4.3, and at redshift z = 0.45. The solid curve denotes a Gaussian radial profile (or off-centered) model with the width that is taken to be r of j = 
0.3r a (~ fOO kpc), where r s is the scale radius of the dark matter NFW profile. The dotted curve is an NFW radial profile model, which is given by the 
concentration parameter of c Q g = 20. 




clusters, and found that the LRG radial distribution can be fitted with an NFW profile with h igh concentration parameter (c ~ 20). Using the 
Subaru weak lensing observations for about 20 X-ray luminous clusters, lOguri et al. I feoich fit an elliptical NFW model to the dark matter 
distribution. They found that, for most clusters, the positional difference between the lensing-inferred mass center and the brightest cluster 
galaxy is well fitted by a Gaussian distribution with width of ~ 100 /i _1 kpc, a scale comparable to the positional uncertainties in the lensing 
analysis. For these clusters, the lensing data is consistent with the DLRGs lying in the center of mass of their host halos. Ho wever, in a few 
cluster s, the DLRGs are clearly offset from the center of the potential with characteristic displacements of ~ 400 /i -1 kpc. Johnston et al] 
1 2007) reached a similar conclusion: most DLRGs are in the centers of their halo; however, a handful are significant displaced. However, the 
results are not yet conclusive due to the limited statistics. In this paper, we employ the following two empirical models for a radial profile of 
DLRGs based on these observational implications: 



p off (r;M) = 



(27r)3/2r o yM) 



exp 



' off 



2rl s (M) 



4ir4 j (coffr/r V ir)(l + c oS r/r vil ) 2 1 



(Gaussian offset model), 
(NFW offset model), 



(12) 



where / = l/[ln(l + c ff) — c ff/(l + Coff)] and the prefactor of each model is determined so as to satisfy the normalization condition 
XT" 1 Aitr 2 dr p «(r) = 1. These profiles are specified by one parameter (r ff or c ff ), but differ in the shape. Note that r v i r is specified as a 
function of halo mass and redshift, and c g differs from the concentration parameter of dark matter profile. 

As a working example, we will employ r s(M) = 0.3r s (M) — 0.3r vlI (M) / c v ir(M) for Gaussian DLRG radial distribution and 
c ff = 20 for NFW DLRG radial distribution as our fiducial models. Here r s is the s cale radius of dark matter NFW profile, and c v i r is the 
concentration parameter. For the following results we will use the simulation results in Duffy et al. (2008) to specify c v i r as a function of halo 
mass a nd redshift. Our fiducial model of r fr = 0.3r s gives r g ~ 100 kpc for halos of 1O 14 M0, consistent with the results in 
(2Ql3>. 

*Jote that the typical offset of the DLRG from the center of the halo potential varies with halo mass. 



Qguri et al 



Fig.[2]shows the Gaussian and NFW models for the radial profile of DLRGs inside a halo of mass M = 10 14 fe _1 M Q and at z = 0.45. 
Note that for our fiducial model, r s ~ 250 ft _1 kpc. These profiles are much more centrally concentrated than a typical dark matter halo of 
this mass scale, represented by an NFW profile with c = 4.3. 

The Fourier transforms of these radial profiles are analytic functions: 



eM-r 2 oS (M)k 2 /2], 



p oS (k;M) = 



f 



sin {Si(r?(l + Coff)) - Si(r,)} + cosr; {Ci(77(l + v )) - Ci(r?)} 



sin(r/c fi 



r/(l + Coff) 



(Gaussian model), 
(NFW model), 



(13) 



where r\ — fcr v i r /c ff , and Si(a;) and Ci(x) are the sine and cosine integral functions. The Fourier transform p s(k; M) 
scales. 



1 on relevant 
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Figure 3. Left panel: Velocity dispersion of DLRGs averaged over DLRG radial profile, as a function of host halo mass at z = 0.45, assuming the Gaussian 
(solid curve) and NFW (dotted) radial profile profiles as in Fig. [2] For comparison the dashed curve shows the velocity dispersion of dark matter. Right panel: 
Redshift dependence of the velocity dispersions of DLRGs and DM for a halo with mass M = 10 14 fe — 1 Mq. 



2.3.3 Velocity dispersion of DLRGs 

In this paper, we simply assume that galaxies at a given radius r have the following ID velocity dispersion that is determined by the mass 
enclosed within the sphere: 

,Ur;M) = ^l. (14) 
Zr 

At virial radius r v i r , the velocity dispersion is determined by virial mass: 

a v (r = r vir ; M) = 472 W s (^j^) ^ (0) ^ d + ^ , (15) 

where we use the defin ition of virial mass given in terms of the overdensity A(z) at redshift z (we use the fitting formulae given in 
Nakamura & Sutol 1997 ). Since an NFW profile has an asymptotic behavior of M(< r) oc r 2 as r -> 0, the velocity dispersion has the 



limit o- v (r; M) — > as r — > 0. 

The velocity dispersion of LRGs is poorly known (see lSkibba et alfcoi ll . for the first attempt). In Appendix lAl we give an alternative 
model of computing the velocity dispersion by assuming an isothermal distribution for the phase space density of DLRGs within a halo, 
where we properly take into account the different radial profiles of DLRGs and dark matter. 

The averaged velocity dispersion of DLRGs within halos of a given mass scale M can be obtained by averaging the velocity dispersion 
(Eq. QH) with the radial profile of DLRGs: 

4 off (M)= / 4itr 2 drp off (r;M)a 2 v (r;M). (16) 
Jo 

This velocity dispersion has an asymptotic limit when all the DLRGs are at the center of each halo: a v . a ff — > when p ff (r) oc 8d (r). Fig. [3] 
plots the velocity dispersion of DLRGs, a v ]0 a (M) , as a function of halo mass M for a fixed redshift (left panel), and as a function of redshift 
z for a fixed halo mass (right), respectively. The velocity dispersion of DLRGs is larger in more massive halos and at higher redshifts. Within 
the same halo, the velocity dispersion of DLRGs is smaller than that of dark matter by 10-20% in the amplitudes, because DLRGs are more 
centrally concentrated. Eq. l !15t implies that the velocity dispersion of both the DLRGs and the dark matter scales as a Vi0 g(M) oc M 1 ' 3 . 

2.4 Redshift-space power spectrum of DLRG and the covariance matrix 

We use our model for the radial distribution of the DLRGs and their velocity distribution to estimate the effect of the offset on the power 
spectrum. We assume the velocity distribution of DLRGs within halos is Gaussian, where the width of the distribution is given by the velocity 
dispersion (Eq. 1141 ): 

1 r 

fl(Ar||;r,M)<*(Ar||) = -= — — exp 

V2na v>of! (r, M) 



dv h (17) 
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Figure 4. The FoG suppression scale (see Eq. 1211 ). a v a ff /aH(z), as a function of the model parameters of DLRG radial profiles, for the Gaussian radial 
profile model (left panel) and the NFW model (right panel). The FoG scale is computed by averaging the velocity dispersion of DLRGs over the halo mass 
function weighted by the DLRG halo occupation distribution (Eq. 1111 ). The different curves are for different redshifts. 



where Ar 1 1 = v^/aH(z) and the prefactor is determined so as to satisfy J^^dAr'^ J?(Ary; M) = 1. 
The Fourier transform of the redshift-space radial profile (see Eq. (6)) can be expressed as 



p s ,off(k±,k\\;M) 



Anr drp «{r; M) exp 



qg, off (r,Af)fcV 
2a 2 H 2 (z) 



(18) 



where k = yj k\ + k 2 and the exponential function above is the Fourier transform of Eq. d!7t . We again note that, exactly speaking, p S)0 s 
also depends on the Fourier transform of the real-space radial profile p s as implied by Eq l[6j, but we use p g ~ 1 at large length scales of 
interest, much larger than the characteristic offset of the DLRG from the halo center. Also note that, for the limit p g — > 8d{t), p s ,off 1 
as o\,, ff — > at r — > 0. 

Hence the redshift-space power spectrum of DLRGs (see Eq. (8)) can be computed for a given cosmological model by inserting Eq. d 1 8b 

into 

1 2 



-Ps.DLRG (k : p) = 



1DLB.G 



rift 

dM —b{M)Nno D (M)p s , oS {k, fx; M) 



At very large length scales (or very small fc's) the redshift-space power spectrum can be approximated as 



f s ,DLRG(fc, pC) ~b 



—2 i 2 2 



i 2 H 2 (z 



Ps,m{k, Ll), 



where a v oS is the velocity dispersion averaged over the halo mass function weighted with the DLRG HOD: 



—2 

°Voff 



brin 



dM^b(M)NuoB(M)al oS (M). 



(19) 



(20) 



(21) 



As we will show below, the approximation l !20t is not accurate at k ^ 0.15ft/Mpc. 

Thus the key quantity characterizing the FoG effect on DLRG power spectrum is the halo-mass averaged velocity dispersion, a% off . 
Tabled gives the values for different redshifts assuming our fiducial model parameters. The DLRG velocity dispersion is also compared with 
that of dark matter within the same halos hosting DLRGs. It can be found that the typical FoG suppression scale, estimated as W V:0 fi/ (aH), 
is of scales of 5 ft _1 Mpc. Therefore, even at large length scales k ~ 0.1 /iMpc" 1 , which is employed in the literature in order to extract 
cosmological information, the FoG effect suppresses the power spectrum amplitudes by a factor of 0.75 (1 — [0.1 x 5] 2 ~ .75) according to 
Eq. d20t . a systematic correction that is much larger than the reported statistical errors in many surveys (see Appendix C in iReid et alfeoid 
for discussion). 

Fig-IDplots how the typical FoG suppression scale changes with changing parameters of the Gaussian and NFW DLRG radial profiles. 
For the Gaussian radial profile, the FoG displacement scale has a maximum scale of a Vt0 g/(aH) ~ 5 /i _1 Mpc around r s ~ r s , as the 
velocity dispersion peaks at the scale radius r s for an NFW profile. On the other hand, for the NFW radial profile, the displacement scale 
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Figure 5. This figure shows how different mass halos contribute to the power spectrum suppression shown in Fig. |4] While the typical DLRG sits in a halo 
with mass ~ 10 13 Mq, the FoG in the more massive halos plays the most important role in suppressing the galaxy power spectrum. The left and right panels 
are the results for the Gaussian and NFW radial profile profiles, respectively, where we assume the fiducial model parameters r Q g = 0.3r s and c g = 20, 
respectively, as in Fig. Ff] 



has a weak dependence on c g, and leads to ~ 5 /i _1 Mpc over all the range of c g. The redshift dependence of the FoG scale is weak and 
changes by only ~ 10%. 

Fig-H]shows that the more massive halos are responsible for most of the FoG suppression. While the typical DLRG sits in halos of mass 
5 x 10 13 r'itf , the FoG arises primarily from more massive halos as these halos have larger velocity dispersions. Here, we again assume 
the fiducial model parameters of the FoG effects as in Fig. [2] The plot shows that halos with masses M ~ 10 14 hT^M® have a dominant 
contribution to the FoG effect at redshifts z — 0.35 and 0.7 for both the Gaussian and NFW radial profiles, while less massive halos become 
more important at higher redshifts. 

The covariance matrix describes statistical uncertainties in me asuring the redshift-spa ce power spectrum from a given survey, and the 
correlations between the power spectra of different wavenumbers. iTakahashi et alj 1(2009) show that the assumption of Gaussian errors is 
valid on scales of interest. In this limit of Gaussian errors, the covariance matrix has a simple form: 

1 



25 5 

Cov[P s ,DLRG(fci, A*a )-Ps, DLRG (fcj, /J-b)] = : ' 



N mode {ki,fi b ) 



Ps,DLRG(fci, fJkb) + 



HDLRG 



(22) 



where ki and fi a are the i-th and a-th bins of wavenumber and cosine angle, respectively, and 5^ and S^ b are the Kronecker delta function: 
6ij = 1 if i = j within the bin width, otherwise 5y = and so on. The Kronecker delta functions impose that the power spectrum of different 
wavenumber bins are independent. The quantity N rnodc (ki, fj, a ) is the number of independent Fourier modes around the bin centered at ki 
and ix a with widths Afc and Afi, which can be resolved for a given survey volume V s : N m "(ki, yu ) = 2-7rfc^ AfcA/iV s /(27r) 3 . Here we 
assume the fundamental Fourier mode is determined by the survey volume as kf — 2n/L, a reasonable approximation for a simple survey 
geometry. 



3 ANGULAR POWER SPECTRUM OF DLRG-GALAXY WEAK LENSING 



Obser vations of DLRG-galaxy lensing measure the radial distribution of DLRGs in the halo. In this subsection, we briefly review lOguri & Takada 
1 201 II) discussion of how the radial distribution of the DLRGs in their halo affects the galaxy-galaxy lensing observables. 

The halos hosting DLRGs distort background galaxy images. By cross-correlating positions of DLRGs on the sky with tangential 
ellipticity component of background galaxy images with respect to the line c onnecting DLRG and background galaxy, we can measure the 
radially averaged mass distribution around a DLRG dMandelbaum et alJ2 006). While this stacking analysis is usually done in real space, we 
will describe the result s in Fourier space as the effect of the DLRG offsets are convolution in real space and multiplication in Fourier space 
l Oguri&Takadalboilb . 

Since we are inte rested in small angular scales, we can use the flat-sky approximation I Limber 1954) and use the halo model in 
Oguri & Takadal 11201 ll) to compute the angular power spectrum of DLRG-galaxy lensing (also see iTakada & Bridl3l2007h : 
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C\ a (l) = C#(J)+C^(J), (23) 

where Cl h n and are t he 1- and 2-halo term spectra defined as follows. For the full-sky expression of the lensing power spectrum, see 
de Putter & Takadal Wd) . The 1-halo term contribution to galaxy-galaxy lensing arises from the mass distribution within one halo that hosts 



DLRGs and gives dominant contribution to the signal on small angular separations: 



C%(1) = f dx-^;SDLRG(z)W GL (x)x~ 2 I dM ^-Nhod(M)— [Mm FW {k; M, z)p oH {k; M)| fc=l/ + m sh , DLRG ] 

n DLRG J «X<«' J O.M p m 

where p m o is the mean mass density today, \ ' s me comoving angular diameter distance (which is given as a function o f redshift via the 
distan ce-redshift relation), W Gh (x) is me lensing efficiency function for a given source galaxy population (see Eq. 19 in Oguri & Takadal 



(24) 



201 1), and d V/dxdft is the volume element in the unit comoving interval and the unit solid angle; d V/dxdQ, = x f° r a nat universe. The 



function Sdlrg(z) is the redshift selection function of DLRGs. For simplicity, we assume a complete selection function: 5*dlrg(z) = 1 
within the redshift range of the survey, and otherwise Sdlrg = 0. The quantity ji-dlrg is me mean angular number density of DLRG in the 
redshift slice: nq^ RG = J dx(d 2 V/dxd£l) J dM(dn/dM)Njior>(M)SY)i,B,G( z )- The term denoted by m s h,DLRG gives the contribution 
arising from a subhalo hosting DLRG, an d we will throughout t his paper assume the subhalo mass ttIsIi.dlrg — 0.32 x 10 12 /i _1 Mq as our 
fiducial value, implied from the results in ljohnston et al.l (l2007h . 

There are two contributions to the 1-halo term in Eq. (1241 : the first term in the bracket describes the contribution of the halo of mass 
M and the second term describes the contribution of a subhalo hosting the DLRG. For the first term, we assume an NFW p rofile char- 
acterizing the dark matter distribution within a halo, and unfw is the Fourier-transform (see Eq. 29 in lOguri & Tak ada 2011). Including 
off-centered DLRGs in the g alaxy-galaxy lensing analysis dilutes the measured lensing signal amplitudes at the small scales l Johnston et alj 
2007 ; Oguri & TakadalEoi 1 ). This off-centering effect on the lensing power spectrum can be included by simply replacing the dark matter 



profile with up s, where p ff is the Fourier- transformed coefficients of the DLRG radial profile (see Eq. 1131 ). For the subhalo contribu- 
tion we simply assume the delta function for the mass profile, a good approximation at the relevant angular scales. In this limit, the power 
spectrum behaves like a white shot noise. 

Similarly, the 2-halo term contribution, which dominates at large scales, is given as 



^ ( °"^y dx S SDLRG(z)H/GL(x)x " 



dM—b{M)Nn OD (M) 
dM 



I 



where P^(fc) is the linear mass power spectrum. This 2-halo term scales with halo bias: if DLRGs are residing on more massive halos or 
equivalently more biased halos, the 2-halo term has greater amplitudes. 

To perform para meter forecasts for planned lensing surveys, we also need to model the lensing power spectrum covariance. Following 
Takada & Jain (2009). we assume Gaussia n errors so that the cov ariance matrix of the lensing power spectrum is given by a product of 



sampling variance and shot noise terms (see lOguri & TakadalboilL for the definition of the covariance matrix). 

Fig. |H shows the angular power spectrum of DLRG-galaxy weak lensing expected when cross-correlating DLRGs in redshift slice 
0.3 < z < 0.4 with background galaxy images that have a typical redshift of z ~ 1 as expected for a Subaru-type imaging survey. The 
top dotted curve shows the power spectrum when all DLRGs are at each halo's center, while the bold solid curves show the spectrum with 
the off-centering effect assuming our fiducial models of the Gaussian (left pane!) and NFW (right panel) radial profiles as in Fig. [2] The 
figure shows that the off-centering effect significantly dilutes the lensing power spectrum amplitudes at angular separations smaller than 
the projected offset scale. The dashed curve shows the contribution from DLRG subhalo assuming m s h,DLRG = 0.32 x 10 12 h' 1 M©. 
The DLRG subhalo contribution becomes significant at I Si 2 x 10 4 , where an angular scale of I ~ 10 4 corresponds to the projected scale 
~ 300 kpc for the DLRG redshift of z = 0.35 for a ACDM model. 

The boxes around the curve show statistical uncertainties in measuring band powers at each multipole bins, expected when measuring 
the DLRG-galaxy lensing for an overlapping area of 2000 square degrees between spectroscopic and imaging surveys. Here we assume the 
depth expected for Subaru HSC survey that probes galaxies at typical redshifts of z s ~ 1 (see Section |4~T| for details). The plot clearly implies 
that combining such imaging and spectroscopic surveys allows us to infer the off-centering effect at a high significance. We will below give 
a more quantitative estimate. 



4 RESULTS 

In this section, we estimate the ability of ongoing and planned surveys for using the DLRG-galaxy weak lensing measurements to correct the 
FoG effect on redshift-space power spectrum of DLRGs. 



4.1 Survey parameters 

To model the DLRG power spectrum available from ongoing and upcoming spectroscopic surveys we assume survey parameters that resemble 
BOSS, Subaru PFS, and Euclid surveys. The expected statistical uncertainties in measuring the DLRG power spectrum in each redshift slice 
depend on the area coverage (or equivalently the comoving volume) and the number density and bias parameter of DLRGs. For the sky 
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Figure 6. Angular power spectra of DLRG-galaxy lensing assuming that spectroscopic and imaging surveys are available for an overlapping area of 2000 
square degrees and the DLRGs are in the redshift range of 0.3 < z < 0.4 (see Section pTTl for details). The bold solid curves are the angular power spectra 
including the off-centering effect for the Gaussian (left panel) and NFW (right) DLRG radial profiles. The thin curves in the left panel are the 1- and 2-halo 
term contributions. For comparison, the dotted curves are the spectra without the off-centering effect. The dashed curves are the contribution from a subhalo 
hosting DLRG, assuming the subhalo mass m s h,DLRG = 0.32 X 10 12 h~ 1 Mq. The boxes around the bold solid curve are the expected la uncertainties in 
measuring band powers of the power spectrum at each I bins, where we assume survey parameters given in Table[T] 



coverage we assume 10,000, 2,000, and 20,000 square degrees for the BOSS, PFS, and Euclid surveys, respectively. Table [T] summarizes the 
survey parameters. For the redshift slices we consider 3 redshift bins with width Az = 0.1 over the range 0.3 < z < 0.6 for the BOSS 
survey, 4 slices with Az = 0.2 over 0.6 < z < 1.4 for the PFS survey, and 10 slices with Az = 0.1 over 0.5 < z < 1.5 for the Euclid 
survey. Thus the BOSS and PFS surveys are complementary to each other in their redshift coverages, while Euclid survey covers a wide 
range of redshifts on its own. Having a wider redshift coverage allows us to trace the reds hift evolution of mass clustering growth over a 
wider range of redshifts and therefore improves cosmological constraints (Taka da et al hood) . In Table Q] we also show the comoving volume 
of each redshift slice. We simply assume the same HOD model given in Eq. {TO over all the redshifts to estima te the number den sity of 
DLRGs and the mean bias. The number densities are somewhat smaller than those actually obtained from data (e.g., White et al.l2011 ) or the 
target number densities for t he survey. This is b ecause we use the HOD model for DLRGs: we assume that we can select one DLRG per each 
halo based on the method of Reid et al. i 20ld) . TableQ]also shows the FoG displacement scale a Vi0 g /(aH) in each redshift slice, which is 
estimated using Eq. 1211 . 

The parameter forecast is also sensitive to the maximum wavenum ber fc max , because complex non-linearities alter the power spectrum 
at high wavenumber. In our analysis, we use Taruva et al. ( 20091 2010h to determine fc max : 



k 2 

'•'max 

6tt 2 



dkP^(k) = C. 



(26) 



We use either C = 0.2 or 0.7, motivated by the fact that the model predictions based on the standard and improved perturbation theory are 
sufficiently accurate in a sense that the predictions well match N-body simulation results to a few % accuracies in the amplitudes up to the 
determined fc max . Table[T]gives the fc max values for C = 0.2 and 0.7, r espectively, in eac h redshift slice. The approach in our analysis is 
to assume that refined model prescriptions based on N-body simulations l Taruva et alfeoioh are used to estimate non-linearities in the real 
space power spectrum and focus on the FoG effects on the DLRG power spectrum. Our goal is to show how the lensing analysis can calibrate 
the FoG effects. 

For our analysis, we assume a lensing survey with either the properties of the planned Subaru HSC survey or the proposed EUCLID 
imaging survey. For the Subaru HSC survey we assume the survey area fi s = 2, 000 square degrees, the mean number density of imaging 



galaxies n™ agmg = 30 arcmin , and the redshift distribution is given by the functional form n g (z) 



£p(— z/zq) where the parameter 



zo is determined so that the mean redshift (z) — 1, For the Euclid survey we simply assume the same survey parameters, except for 
the survey area of fi 8 = 20, 000 square degrees. When computing the DLRG-galaxy lensing we use background galaxies at redshifts 
z > ZDLRG.i + Az/2 + 0.05, where zdlrg.i + Az/2 is the upper bound on redshifts of DLRGs in the i-th redshift slice of a given 
spectroscopic survey. That is, we include galaxies at redshifts higher than the redshifts of any DLRGs in a given r edshift slice by 5z = .05. 
We assume that such background galaxies can be selected based on their available photometric redshift estimates ( Nishizawa et alj201fj) . In 
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[(Gpc/ft) 3 ] 


[ft/Mpc] 


[h/Mpc] 


[(/i/Mpc) 3 ] 


b 


[10 14 Af Q /h] 


DM 


DLRG 


DLRG 






















(Gauss) 


(NFW) 


BOSS 


0.3- 


0.4 


0.76 


0.10 


0.17 


1.25 x 10~ 4 


1.80 


0.63 


5.35 


4.42 


5.06 


(10,000 deg 2 ) 


0.4- 


0.5 


1.14 


0.10 


0.18 


1.15 x 10~ 4 


1.85 


0.58 


5.41 


4.48 


5.10 




0.5- 


0.6 


1.53 


0.11 


0.19 


1.05 x 10~ 4 


1.89 


0.54 


5.43 


4.50 


5.11 


PFS 


0.6- 


0.8 


0.84 


0.11 


0.20 


0.91 x 10~ 4 


1.95 


0.48 


5.42 


4.49 


5.07 


(2,000deg 2 ) 


0.8- 


1.0 


1.13 


0.12 


0.22 


0.72 x 10~ 4 


2.02 


0.42 


5.33 


4.43 


4.97 




1.0- 


1.2 


1.37 


0.13 


0.23 


0.57 x 10~ 4 


2.08 


0.36 


5.20 


4.32 


4.81 




1.2- 


1.4 


1.56 


0.14 


0.25 


0.44 x 10~ 4 


2.14 


0.31 


5.03 


4.19 


4.64 


Euclid 


0.5 - 


0.6 


3.06 


0.11 


0.19 


1.05 x 10~ 4 


1.89 


0.54 


5.43 


4.50 


5.11 


(20,000deg 2 ) 


0.6- 


0.7 


3.84 


0.11 


0.20 


0.95 x 10~ 4 


1.93 


0.50 


5.56 


4.60 


5.09 




0.7- 


0.8 


4.60 


0.12 


0.20 


0.86 x 10~ 4 


1.96 


0.46 


5.53 


4.59 


5.05 




0.8- 


0.9 


5.32 


0.12 


0.21 


0.77 x 10~ 4 


2.00 


0.43 


5.49 


4.56 


5.00 




0.9- 


1.0 


5.98 


0.13 


0.22 


0.69 x 10- 4 


2.04 


0.40 


5.43 


4.51 


4.93 




1.0- 


1.1 


6.58 


0.13 


0.23 


0.61 x 10" 4 


2.07 


0.37 


5.36 


4.45 


4.85 




1.1 - 


1.2 


7.11 


0.14 


0.24 


0.54 x lO" 4 


2.10 


0.34 


5.28 


4.39 


4.77 




1.2- 


1.3 


7.58 


0.14 


0.25 


0.47 x 10~ 4 


2.13 


0.32 


5.20 


4.33 


4.68 




1.3- 


1.4 


8.00 


0.15 


0.26 


0.41 x 10~ 4 


2.16 


0.30 


5.11 


4.26 


4.59 




1.4- 


1.5 


8.35 


0.15 


0.26 


0.36 x 10" 4 


2.18 


0.28 


5.01 


4.18 


4.50 



Table 1. Survey parameters considered in this paper, which are chosen to resemble the ongoing BOSS survey, and the planned Subaru PFS survey and Euclid 
survey. The range of redshift z, the survey volume of each redshift slice Vs, the mean number density of DLRGs Sdlrg. the effective linear bias b, and the 
mean halo mass are given. The quantity fc mac ,C=0.2 or ^mac,C=0.7 is me maximum wavenumber of each redshift slice up to which the power spectrum 
information is included when studying parameter forecasts. The maximum wavenumbers are chosen by setting C = 0.2 or 0.7 in Eq. {26), which is motivated 
by the fact that there are accurate model predictions of halo clustering available up to such scales. The last three columns denote the FoG suppression scale, 
~®v,os/ a H(z) [Mpc/h], for dark matter and for DLRGs having the Gaussian and NFW radial profiles. 



addition we assume the rms intrinsic ellipticities a e = 0.22 per component , which determines th e intrinsic ellipticity noise contamination to 
the error covariance matrix of the DLRG-galaxy lensing power spectrum ( Oguri & Takada 201 1 ). 



4.2 Weak lensing information on the off -centered DLRGs 

In this subsection, we estimate the ability of future lensing survey to determine the radial distribution of DLRGs in the halos. 

As shown in Fig. [6] the offset of the DLRGs from the centers of their halos dilutes the lensing signals at small angular scales. Thus if 
we assume that the dark matter halos are well-described by NFW halos on the ~100 kpc scale, the lensing observations can be used to infer 
the DLRG radial distribution. 

There are two approaches that can be considered. The first one is a method fully based on the halo model. That is, by fitting the measured 
lensing profile to the halo model prediction (see Eq. 11241 ). one can constrain parameters including the parameters of either Gaussian or NFW 
DLRG radial profile model: 



{r off (Gauss) orc off (NFW),Civ,/?,(Ti ogA /,m shiDLR G} , 



(27) 



where r a g or c D s is the off-centering parameter (see Eq. 11121 ). eri og m is the parameter of DLRG HOD (see Eq. 1111 ). and the mean mass scale 
of subhalo hosting each DLRG (see Eq. 1241 ). The parameters cjv and /3 are introduced to model the scaling relation of halo concentration 
with halo mass, c v ir(Af v ir) = Cjv(M v i r /10 12 ft. — 1 Mq)~^ , where cjv = 7.85 and /3 = 0.081 are adopted for the fiducial values. Note that 
the HOD models are given by the two parameters <ri og m and M m i u , but one of the two is determined by the observed mean number density 
of DLRGs. In this parameter estimation we assume that the halo mass function, halo bias and NFW profile or more generally halo mass 
profile are well calibrated based on a suite of simulations. We also assume that, although the lensing strength depends on cosmology, the 
background cosmology is well constrained from other cosmological observables such as CMB and the galaxy redshift survey itself (e.g., via 
the BAO experiment). Hence we include only the 1-halo term power spectrum of DLRG-galaxy lensing to perform the parameter forecast 
based on the Fisher matrix formalism. 

We can use ex ternal information as priors on the model paramet ers. For the concentration-mass scaling relation, simulation-based studies 
l Duffy et al.l2008h and/or cluster lensi ng studies Jokabe et al.l20ld) can be used as the priors. For the HOD parameter, ai og m, the clustering 
analysis of DLRGs i White et al]l2011 ) can be used. We employ the priors cr(lncjv) = a (In/?) = a (In <ri og m) = 0.2 in our Fisher analysis. 

Table |2] gives the marginalized errors on the parameter, r g or c ff, for the Gaussian or NFW radial profile model, respectively, in each 
redshift slice. Here we assume the DLRG-galaxy lensing measurements for an area of 2,000 or 20,000 square degrees for the Subaru HSC or 
Euclid survey, respectively, as discussed in § 14.11 Note that we include the lensing information up to i max = 5 x 10 4 in our Fisher analysis, 



© 0000 RAS, MNRAS 000, 000-000 



Using galaxy-galaxy weak lensing measurements to correct the Finger-of-God 13 



Off-centering Parameter Determination from DLRG-Galaxy Lensing 







Gaussian offset 


NFW offset 


Survey 


Redshift range 


o-(ln r off ) 


""(fncr^off) 


cr(lnc ofI ) 


cr(ln cr„, ff) 


BOSS 


0.3-0.4 


0.07 


0.019 


0.31 


0.016 


(10,000 deg 2 ) 


0.4-0.5 


0.08 


0.022 


0.33 


0.018 




0.5 - 0.6 


0.10 


0.027 


0.34 


0.019 


PFS 


0.6-0.8 


0.12 


0.032 


0.34 


0.020 


(2,000 deg 2 ) 


0.8 - 1.0 


0.24 


0.066 


0.36 


0.022 




1.0- 1.2 


0.49 


0.15 


0.39 


0.024 




1.2-1.4 


0.85 


0.33 


0.46 


0.030 


Euclid 


0.5 - 0.6 


U.UJO 


U.UIO 




n ni 9 


(20,000 deg 2 ) 


0.6-0.7 


0.087 


0.024 


0.24 


0.014 




0.7-0.8 


0.13 


0.036 


0.27 


0.015 




0.8 - 0.9 


0.20 


0.055 


0.29 


0.017 




0.9-1.0 


0.31 


0.086 


0.31 


0.019 




1.0-1.1 


0.47 


0.14 


0.33 


0.020 




1.1 - 1.2 


0.71 


0.24 


0.35 


0.022 




1.2-1.3 


1.03 


0.51 


0.37 


0.023 




1.3-1.4 


1.44 


0.53 


0.39 


0.025 




1.4-1.5 


1.87 


0.54 


0.41 


0.027 



Table 2. Marginalized uncertainties in the DLRG off-centering parameters in each redshift slice, which are expected to obtain from the DLRG-galaxy lensing 
measurements (see around Eq. 1271 for details of our Fisher analysis). For the BOSS and PFS surveys, we assume that the Subaru HSC-type imaging survey is 
overlapped with the spectroscopic surveys for an area of 2000 square degrees. For the Euclid survey we assume the overlapping area of 20000 square degrees 
for the joint imaging and spectroscopic surveys. As in Fig.|2]we consider the Gaussian or NFW DLRG radial profiles as a working example, and the the error 
bars show the relative errors: e.g., <r(lnr j}) = &(r f{)/r ff , where the denominator is the fiducial value. The relative error on the FoG suppression scale, 
denoted by cr(ln a v g ), is estimated by propagating the error of r g or c g using Eq. 12 It . 



but we find that the results in Table[2]are not so sensitive to the choice of the maximum multipole as such small-scale signals are dominated 
by the lensing contribution of DLRG's subhalo. Most of the information on the characteristic offsets comes from the power spectrum at 
I 10 4 . Table|2]also shows the expected error on the FoG suppression scale, which is obtained by propagating the error of the characteristic 
offset based on Eq. 12111 . 

The method above is a model-dependent method. Since the Gaussian or NFW radial profile models may differ from the genuine radial 
profile, the estimated offset parameters may have systematic biases. Moreover the radial profile may change with halo masses. Hence an 
alternative approach is based on a metho d that the measured l e nsing profile is compared with an single NFW profile, including the off- 
centering effect. Previous measurements ( Johnston et al. 2007 ; Okabe et al. 2010h mostly used this method. However, we find that this 
method does not work well if we use the full sample of DLRGs for the DLRG-galaxy lensing measurements. If the lensing profile in Fig.[6]is 
fitted to a single NFW profile, the extracted halo mass is found to be ~ 0.5 x 10 14 h' 1 M Q , i.e. the mean halo mass of halos hosting DLRGs 
as listed in TableQ] This mass scale is smaller than the mass scale of halos, 10 14 /i _1 M Q , that give a dominant contribution to the FoG effect 
as shown in Fig. [5] In turn the characteristic offset inferred from such a single NFW profile-fitting is not that accurate. Since the FoG effects 
are larger for the more massive halos, it would be more optimal to focus on studies of the DLRG-galaxy lensing for halos with a narrower 
mass range around lO 14 /i -1 M0, where such a halo catalog can be constructed based on optical richness available form multi-color imaging 
survey itself, X-ray and/or Sunyaev-Zel'dovich (SZ) data. 

In fact we are now working on the SDSS DR7 LRG catalog in order to explore the feasibility of the method above with real data. In this 
on-going study, we have focused on the regions including multiple LRGs in the small spatial region, which are likely to reside in the same 
halo with masses more massive than other majority of halos hosting a single LRG inside. We have then studied the galaxy-galaxy lensing 
signal measured via cross-correlation of the LRG regions with background galaxy shapes, where the background galaxies are taken from the 
photometric SDSS galaxy catalog based on their photometric redshift estimates. We have so far found, preliminarily though, that the lensing 
signals at small scales do change with different centers; the brightest LRG position, the faintest LRG or the mean of their positions. Then we 
are now trying to constrain the off-centered profile of LRGs from the varied lensing signals including marginalizations over other parameters 
such as the mean halo mass, the halo profile parameters and the sub-halo mass scale. This is in working progress, and will be presented 
elsewhere (Hikage et al. in prep.). 

In this paper, we will just assume that the first approach stated above, i.e. the halo model based fitting method, is feasible for the 
following analysis. 
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Figure 7. Left panel: Redshift-space power spectrum of DLRGs, P^dlrgC^ ' k±)> m tne two-dimensional wavenumber space (fcy , k±). We use Eq. (8) 
to compute the power spectrum assuming the Gaussian DLRG radial profile in Fig. [2] and assume the parameters expected for BOSS survey in Table [T] 
The spectrum shown is the spectra averaged over the redshift range 0.3 ^ z ^ 0.6. The contours denote the power spectrum amplitudes of 10 5 ,3.2 X 
10 4 , 10 4 , 3.2 X 10 3 , 10 3 , 3.2 X 10 2 , 10 2 (h/Mpc) 3 , respectively. Right panel: The monopole power spectrum of DLRGs obtained by averaging the redshift- 
space power spectrum over the angle fi(= fc||/fc) for a fixed k(= yj^± + k 2 ). For comparison, the non-linear matter power spectrum P^ L (fc) and the 

real-space power spectrum of DLRGs (Eq. 0) are plotted. The error bars show the la statistical errors in measuring the power spectrum at each k bin for 
BOSS survey. 



4.3 The FoG effect on DLRG power spectrum 

In this subsection, we explore how the off-centered DLRGs suppress the power spectrum through their FoG effects. The left panel of Fig. [7] 
shows the DLRG power spectrum P s ,dlrg (k± ,k») (Eq. (8)) in the two-dimensional space of (k± ,ku) assuming the Gaussian DLRG radial 
profile model. The plotted power spectrum is the spectra averaged over a range of redshifts, 0.3 ^ z ^ 0.6, covered by the BOSS survey. 
Apparent anisotropic features can be found: at large length scales (k 0.15/i/Mpc), the Kaiser effect enhances the power along the line-of- 
sight direction, while the FoG effect suppresses the power more significantly with increasing wavenumbers and thus squashes the contours 
along fen. The right panel shows the monopole power spectrum which is obtained by averaging the redshift-space power spectrum over the 



angle /i (= fey/fc) for a fixed fc(= y k 2 + k\). For comparison we also show the real-space DLRG power spectrum Pdlrg (&) (Eq. J2]) 
and the input non-linear matter power spectrum P„ L (fc), where the latter is computed using the fitting formula bv lSmith et al.l ( 12003b . It is 



clear that the redshift-space distortion effect causes a scale-dependent modification in the amplitude as well as shape of monopole power 
spectrum. 

To be more quantitative, Fig. [8] shows the ratio between the monopole spectra of redshift- and real-space DLRG power spectra as a 
function of k. Since the linear Kaiser redshift distortion causes an overall offset in the monopole power spectrum amplitudes, the scale- 
dependent effect in the figure is solely due to the FoG effect. The data points show the DLRG power spectrum assuming the Gaussian radial 
profile model. The error bars are the expected 1 a statistical uncertainties in measuring band powers of the power spectrum at each k bin 
for the BOSS survey, which include the sampling variance and shot noise contamination. The plot implies that the FoG effect can be very 
significant even at very large length scales, k ~ 0.1 /iMpc -1 , relevant for BAO scales, and the FoG suppression is more significant at larger 
fe's. The FoG effect suppresses the power spectrum amplitudes by 10 (30%) around k = 0.1(0.2) /iMpc -1 . 

The following empirical models of the FoG effect are often assumed in the literature: 

Gaussian : exp[— (k)j,a Vi0 s /aH(z)) 2 ], 

Lorentzian : _ 1 , (28) 

1 + {k^a VtOS /aH{z)Y 

which are called the Gaussian and Lorentzian FoG models. These two models are given by a single parameter a V:0 S an d have the same form 
given as 1 — (kfia Vi0 f[ /aH) 2 at the small k limit. Note that we use the same notation W Vl0 e as that of the radial profile model for notational 
simplicity, but keep in mind that a v ,as in the equation above is a free model parameter. 

The dashed and dotted curves in Fig. UJ show the results for these FoG models, which can be compared with the FoG effect due to off- 
centered DLRGs we have developed in this paper. Here the model parameter <?„ >0 g is taken from Eq. d2]j assuming the input Gaussian radial 
profile model. At scales k £ 0.15/iMpc - , all the results well agree with each other because the FoG effect can be well approximated as 
1 — (fc/i<T„, ott - / aH) 2 at such small fc's. At the larger k's the model differences become significant: the Gaussian approximation underestimates 
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Figure 8. Left panel: The relative difference between the monopole redshift-space power spectra of DLRGs with and without redshift distortion effect (the 
linear Kaiser effect plus the FoG effect). The data at each k bins shows the halo model predictions we have developed in this paper, and the error bar at each k- 
bin shows the lcr statistical uncertainties of the band power expected from the BOSS survey (see Tablefj}. For comparison, the dotted and dashed curves show 
the model predictions assuming a Gaussian or Lorentzian FoG form (see Eq. 1281 ). which is specified with a single parameter cr off (Eq. 1211 ). The full halo 
model prediction and the approximations show differences greater than the statistical error, more significantly at higher fe's. Right panel: The difference ratio of 
the power spectra between the full FoG effect and the approximated FoG effects (Gaussian or Lorentzian FoG models): qlBjq/^» DLR.G — ^' Also snown 
is how the redshift-space power spectrum changes with changing the dark energy equation of state parameter wq and the neutrino mass f v {= Cj/o/f^mo): 
here we consider variations to too = —0.9 and /„ = 0.01 (corresponding to the total neutrino mass m^tot = 0.104 eV) from the fiducial values wo = —1 
and /„ = 0, respectively. 



the power, while the Lorentzian overestimates the power. The differences arise because the FoG effect in our model arises after the integration 
of halo mass function and the DLRG radial profile (see Eqs. 1181 and 1191 ), and therefore the resulting FoG effect does not exactly follow a 
Gaussian form, even if we assume that the velocity probability distribution of DLRGs is Gaussian for each halos. 

The right panel of Fig.[8]more explicitly shows the relative differences between the FoG effect and the FoG approximations (Eq. 1281 ). 
The plot shows that either Gaussian or Lorentzian FoG model cannot be accurate enough compared to the accuracies of upcoming galaxy 
surveys. The plot also shows how the monopole power spectrum changes with changing the dark energy equation of state parameter to 
wo — —0.9 from wq — —1 or the neutrino mass to /„ = 0.01(m„, to t = 0.104 eV) from f v = 0. These parameters are both sensitive to 
the power spectrum amplitudes and therefore most affected by the FoG uncertainty as we will study below more extensively. The plot shows 
that the upcoming surveys have a much higher statistical precision than the effects due to these parameter changes, and also shows that a 
correction/calibration of the FoG effect is very important in order not to have a biased estimate on these parameters. Finally, for the FoG 
effect due to the NFW radial profile model, we have also found similar results to the results in Fig. [8] 

The amount of the FoG suppression effect shown in Fig. [8] obviously depends on the fiducial parameters of our FoG model. One of 
the model uncertainties is the velocity distribution of DLRGs within a halo. If we use the model given in AppendixfA] the FoG suppression 
effect is roughly half of the amplitude shown in Fig. [8] A full analysis of both simulations and observations will be essential to quantify the 
amplitude of this suppression. In the following analysis, for simplicity we will continue to assume the FoG effect computed using our fiducial 
Gaussian and NFW radial profile models. 



4.4 The Impact of FoG Effect on Cosmological Parameter Estimations 

How does an imperfect modeling of the FoG effect affect cosmological parameter estimation? In particular how can adding the DLRG- 
galaxy lensing information helps correct for the FoG effect on the redshift-space DLRG power spectrum? In this subsection, we address 
these questions. 

We perform a Fisher analysis to estimate expected accuracies of parameters for upcoming galaxy surveys, and possible biases due to an 
imperfect modeling of the FoG effect. As the observable we use the redshift-space DLRG power spectrum. In the galaxy clustering analysis 
one needs to assume a reference cosmological model to infer the spatial position of each galaxy from the observed redshift and angular 
position. However, the assumed cosmology generally differs from the underlying true cosmology, which cau ses cosmological distortion 
effect on the observed power spectrum. Hence the observed redshift-space power spectrum can be given as ( Alcock & Paczvnski 19791 ; 
Seo & Eisensteinll2003h : 
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l I <" 1 < - 1 i -P^( 2: )fld H{z) 

Ah Da(z)2 H (zu Pa ' DLRGik "> k±) > (29) 

where fcff d = {D A /D A ,&d)k h k^ = (H &A /H)kx, and D A (z) and H(z) are the angular diameter distance and the Hubble expansion rate 
at redshift z. The quantities with subscript "fid" denote the quantities of the assumed reference cosmology. We use Eq. ( 1191 to compute the 
redshift-space power spectrum. 

Using the power spectrum covariance (Eq. 11221 ). the Fisher information matrix of the DLRG power spectrum measurement is computed 



v ,i ,* m « ainP° b D s LRG (fc,M;^)9inP° b D s LRG(fc,M;^) 



^ LRG d/M I k z dk 



n D LRG(2i)-P s °,DLRG( fc . W Z i) 



(30) 



where the summation is over redshift slices, V Z( is the comoving volume of the i-th redshift slice, and p a denotes a set of model parameters. 
The redshift ranges of each slice are given in Table Q] For each redshift slice, the minimum wavelength kmin is set to be 27r/V z 1 / 3 and the 
maximum wavelength fc max are chosen using the criterion given by Eq. (1261 as listed in TableQ] As we discuss around Eq. ( 1261 . the impact of 
the FoG effec t becomes more si gnificant with including the power spectrum information up to the higher fc max . We use a publicly available 
code CAMB jLewis etal]~2 000) to compute the input linear-mass power spectrum. 



The galaxy power spectrum alone cannot determine all the cosmological parameters due to severe parameter degeneracies I Taka da et al 



2006). Therefore, we combine CMB information with the galaxy power spectrum, which helps efficiently in breaking the parameter degen- 



eracies. We use the CMB information expected from the Planck experiments: the temperature spectrum, _E-mode polarization spectrum and 
the cross-spectrum over multipole range of 2 ^ I ^ 1500. We compute the CMB Fisher matrix F^ anck assuming the instrument noise and 
beam size of each frequency channel in the Planck websitj^i The Fisher matrix for a joint experiment of the CMB and galaxy power spectra 
is simply estimated as F a p — i 7, ( ^) anck + F^ nG . The marginalized error of the a-th parameter, a(p a ), is estimated as (F _1 )q{» 2 , where 
F - is the inverse matrix of the Fisher matrix F. 

The parameter estimation is sensitive to a set of model parameters as well as a choice of the fiducial model. We include a fairly broad 
range of model parameters: 

p a = {fl h h 2 ,Q cdln h 2 , Qk, Hoe, r, A, n s ,l/n(zi),b(zi), r a{zi) or c ff(zi) or a v , ti(zi), /„ and/or w or/ 9 (z = z;)}. (31) 



The fi rst 7 para meters are cosmo logical parameters, and the fiducial values are given at the end of Section [TJ Following Seo & Eisensteinl 



1 2003 ) (also see Saito et aklboil ). we include \jn(zi) and b(z,i) as free parameters in order to model uncertainties due to the residual shot 



noise and galaxy bias in each redshift slice. The FoG effect fully computed based on the halo model is parametrized by the off-centering 
parameter in each redshift (Eq. 11121 ): r g(zi) for the Gaussian radial profile model, or c a g (zi) for the NFW model. When we approximate 
the FoG effect with the Gaussian or Lorentzian form (Eq. 11281 ). we instead use ~a Vt0 ft{zi) (Eq. 1191 ) as a parameter. The fiducial value of 
&v,ott(zi) is computed based on the halo model (eq. 1211 ). 

In Eq. ( 1311 . we further include additional parameter(s): the dark energy equation-of-state parameter too (— 0) and/or the neutrino energy- 
density fraction /„ = Q„/fi m (= O-Olf^l or the growth rate of each redshift slice f g (Zj) = d In Dj In a\ z . The values in the parenthesis 
denote the fiducial values, and the fiducial values of f g in each redshift slice are taken from those of our fiducial ACDM model. These 
parameters are all sensitive to the small-scale amplitudes of galaxy power spectrum, and therefore degenerate with the FoG suppression 
effect. We will pay special attention to an issue of how a knowledge of the FoG effect helps constrain these parameters and minimize a 
possible bias in the parameter caused by the FoG uncertainty. The dimension of our Fisher matrix is at most 19 x 19, 23 x 23, and 47 x 47 
for BOSS, PFS and Euclid surveys, respectively. 

In the parameter forecast we will consider the following four cases for the treatment of the FoG effect: 

(i) The redshift-space power spectrum of DLRGs, measured from a hypothetical survey (BOSS, PFS or Euclid), is compared to the 
model power spectrum, which is given by the non-linear redshift-space matter power spectrum (Eq. J9|) multiplied with either Gaussian 
or Lorentzian FoG model (see Eq. ||28l ). In this fitting, the parameter a Vt0 g of the FoG model is treated as a free parameter (we will hereafter 
call "w/o offset"). 

(ii) The similar approach to the case (i), but we employ the external information on a Vj0 s from the DLRG-galaxy lensing measurements 
assuming combined imaging and spectroscopic surveys: the Subaru HSC survey combined with either BOSS or PFS survey or the combined 
Euclid imaging and spectroscopic surveys (we will call "with offset"). To be more precise we use the statistical error on <r v>0 g in Table|2]as 
the prior of the Fisher analysis. 

(iii) This is fully based on the halo model approach: the FoG effect is computed based on the halo model (Eq. 1191 ) assuming either Gaussian 
or NFW radial profile model. We use the lensing-derived constraints on r g (Gaussian) or c a (NFW) listed in Table [2] as the prior of the 
Fisher analysis (we will call "FoG shape known"). 



12 http://www.sciops.esa. int/index.php?project=PLANCK 

13 The fiducial value f v = 0.01 corresponds to the total neutrino mass m^tot = 0.104eV for our fiducial ACDM model. This value is close to the lower 
limit if the neutrinos obey the inverted mass hierarchy. 
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(iv) This is the extreme case that we neglect the FoG effect in the parameter estimation by setting cr„ i0 ff to be = (we will call "FoG 
neglected").. 



Based on the Fisher matrix formalism, we estimate the bias due to an imperfect modeling of the FoG effect, Sp a (= p a , cs t — Pa, true 
the difference between the estimated and input values of the a-th parameter p a . The amount of the bias is estimated as 



Sp a = 2J(F l ) a pb/3, 
where 

6 " s 



). 

(32) 



1 pk„ 

1 J fcmi 



Ps, DLRG (fc , M Zj) 

(fc,M; z») 



papp 
s.DLRG 



dpf) 



1 + 



nDLRGft,DLRG(^ A*5 Zi 



(33) 



Here P s ,dlrg denotes the true power spectrum where the FoG effect is computed based on the halo model. On the other hand, P^^q 
is the approximated spectrum using Gaussian or Lorentzian FoG model. The total error including both statistical and systematic errors is 
estimated as 



[Ap a ] 2 = [cr s tat(p a )] 2 + (Sp a ) 2 



(34) 



We estimate the impact of the FoG effect on the measurement of the growth rate in each redshift slice. Figs. [9] and \10\ compare the 
marginalized, fractional errors of f g (zi) for the cases (i) and (ii), the cases with and without the FoG effect correction using the DLRG- 
galaxy lensing information of the DLRG radial profile (see TableO. The upper panels show the results for the Gaussian radial profile model, 
while the lower panels show the results for the NFW radial profile model. The left- and right-side panels are different in the maximum 
wavenumber fc max , which is determined by C = 0.2 (left panels) and C = 0.7 (right) in Eq. l !26t roughly corresponding to fc ma x — 0.1 
and 0.2 /iMpc -1 respectively (see Table|2j- Note that the plotted errors include the statistical and systematic errors (Eq. J34]). It is found 
that the lensing information significantly improves the constraint on f g over a wide range of redshifts z: the improvement is up to a factor 
of 2 for C = 0.2 (the lower fcmax), but less significant at higher z. This result can be understood as follows. The DLRG-galaxy lensing 
measurements are more accurate for lower z DLRGs because of higher number densities of background galaxies, which reduces the shot 
noise contamination, (see Table[2}. The error of f g improves less at higher fc where the systematic error due to the FoG model inaccuracy is 
dominated over the statistical error. The Gaussian FoG model gives a better improvement than the Lorentzian model because we assume the 
Gaussian velocity distribution of DLRGs within halos. Finally, comparing the upper- and lower-side panels shows that the NFW radial profile 
model gives a better performance of the lensing FoG correction, as implied by the accuracies of the DLRG velocity dispersion reconstruction 
in Table [2] As a result these spectroscopic surveys allow for constraints on the growth rate to 5% precision or even better at each redshift 
when the lensing information is combined. 

Tables [5] and [4] study the impact of FoG effect on the parameter f v and wo, respectively. The tables show the marginalized error on 
each parameter and the amount of bias in the best-fit parameter if the FoG effect is corrected imperfectly or ignored. Note that in this 
analysis, either wo or /„ is fixed to the fiducial value: for the table of wo, f v is fixed to / „ = and vice versa . Although the dark energy 
equation-of-state parameter wo can be more robustly constrained by the BAO peak location dEisenstein et al ■l2005h , adding the amplitude and 
shape information o f power spectrum can significantly improve the error because dark energy is sensitive to the power spectrum amplitude 
via the growth rate dSaito et al.ll201 If) . First of all, the table shows that the FoG effect is very significant even at such large length scales 
(fcmax — 0.1 /iMpc -1 for C = 0.2). For example, the column "FoG neglected" indicates that the parameters are significantly biased if the 
FoG effect ignored. The bias can be larger than the statistical error by more than a factor of 5. That is, if the FoG effect is ignored, a non-zero 
neutrino mass or Wo — 1 may incorrectly be claimed. 

The table also shows that the parameter accuracy of /„ or wo can be improved by correcting for the FoG effect based on the lensing 
information of DLRGs radial profile. For C = 0.2 (fcmax — 0.1 /iMpc -1 ), comparing the columns labeled as "w/o offset" and "with offset" 
shows that the lensing FoG correction improves the marginalized error of /„ or wo by about 5-25% for BOSS, PFS or Euclid. However, 
as the power spectrum information up to the higher fcmax is included, the bias due to the imperfect modeling of FoG effect becomes more 
significant, as shown in the rows denoted by "C = 0.7". Thus even if the statistical precision is apparently improved in such high-fc regime, 
a more accurate modeling of the FoG effect is needed in order not to have any significant parameter bias. For BOSS or PFS surveys, the 
amount of the bias on /„ or wo is still smaller than the (marginalized) statistical error, while the bias becomes significant for Euclid survey. 
Finally we remark on the results that directions in the parameter bias are opposite for the results with and without the lensing priors on 
the off-centering parameters. For example, for a Gaussian FoG model case, we have found that, if including the lensing prior on the FoG 
suppression scale (i.e. the case with "with offset"), the model prediction tends to overestimate the FoG suppression as implied in Fig. [8] and 
thus prefers a smaller /„, i.e. <5/„ < 0. On the other hand, if the FoG suppression scale is treated as a free parameter (i.e. the case "w/o 
offset"), the fitting tends to prefer a smaller FoG suppression, and in turn prefer a larger neutrino mass (larger /„); 5f v > 0. 

The Gaussian and Lorentzian approximations are not sufficiently accurate to fully describe the FoG shape at fc ^ 0.15/t/Mpc. Fig. 1111 
shows the systematic impact on the marginalized errors of (fu,u>o) for BOSS and PFS surveys with different fcmax values. In this Fisher 
analysis both /„ and too are treated as free parameters in the model fitting (the growth rate f g {zi) is fixed to the fiducial value). The error 
ellipse in each panel shows that the parameter u>o and /„ are correlated with each other in the measured DLRG power spectrum. The figure 
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Figure 9. The marginalized error on the growth rate / 9 (= din D /din a) at each redshift slices of BOSS and PFS surveys, expected when either of the two 
surveys is combined with the Planck CMB information. Note that the error shown here includes both the statistical and systematic error contributions (see 
Eq. 1341 ). The left- and right-side panels differ in the maximum wavenumber fe ma x that is determined by the criterion of either C = 0.2 or 0.7, respectively 
(see Eqs. 1301 and 1261 and Table[TJ. For the upper- and lower-side panels, we assume Gaussian and NFW radial profiles of DLRGs respectively, as in Fig.0 
The long- and short-dashed curves show the results when the FoG effect is modeled by either Gaussian or Lorentzian form and the model parameter a v D ff 
(Eq. 1281 ) is treated as a free parameter in the model fitting. The solid and dotted curves show the results including the DLRG-galaxy weak lensing information 
on a v jf in Table[2] An inaccuracy of the approximated FoG form biases the growth rate estimation more significantly when including up to the higher k's, 
as implied in the right-side panels. 



shows that an imperfect modeling of the FoG effect may bias the parameters. However, the amount of the bias is still at l-a confidence level 
in the two-dimensional space for these surveys. 

Fig- El shows the forecasts for the combined experiments of Planck CMB combined with BOSS, PFS and HSC (upper panels) and 
Planck CMB with Euclid (lower), respectively. The figure shows that the error ellipses are further shrunk by having the power spectrum 
information over a wider range of redshifts. As a result, the marginalized errors of cr(f v ) ~ 0.011(0.009) [or equivalently o"(m„ it ot) — 
0.11(0.093) eV] and a(w ) ~ 0.063(0.051) can be achieved at C = 0.2(0.7) for the BOSS+PFS+HSC survey. In Euclid survey, the error 
decreases to cr(/„) ~ 0.0035(0.0028) [cr(m Vl tot) = 0.035(0.028) eV] and a(w ) ^ 0.029(0.024) at C = 0.2(0.7). In particular the 
Euclid survey allows for a detection of the non-zero neutrino masses at more than lcr significance, because the lower mass bound for the 
normal or inverted mass hierarchy is about 0.05 or 0.1 eV, respectively. Thus such a high-precision determination of these parameters is 
potentially feasible for these surveys ifwe can use a sufficiently accurate model of the FoG in order to include the power spectrum amplitude 
information up to such higher k's. 
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Figure 10. Similar to Fig[9] but for the Euclid survey. 



Marginalized Error and Bias of Neutrino Mass (/„) 
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Table 3. Marginalized error and bias of /„ for either BOSS, PFS, or Euclid survey combined with the Planck CMB information. Here we assume the FoG 
effect due to our fiducial Gaussian DLRG radial profile model, as in Fig.|9] and assume f v = 0.01(T7i„,tot = 0.104 eV) as the fiducial values. Each column 
shows the statistical error (cr(fu)) and the amount of bias = f u> est — /^,truc) due to the difference of the assumed Gaussian (Gauss) or Lorentzian 
(Lorentz) FoG model from the input FoG effect that is computed based on our halo model. The columns labeled as "with offset" and "w/o offset" show the 
results with and without the lensing information being used to correct the FoG effect. The percentage numbers in parentheses in the cr(/„) column denote an 
improvement in the errors due to the lensing information. The column "FoG shape known" shows the statistical error when the same lensing information is 
added but assuming the shape of the FoG effect is known based on the halo model (see text for details). The column "FoG neglected" shows the parameter bias 
where the FoG effect is completely ignored (i.e. ~n v Q ff = is set). These results show a significant bias, more than 100% bias compared to the input value 
/„ = 0.01. 
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Table 4. Similar to Table [5] but for the dark energy equation-of-state parameter wq. Note that for these results the neutrino mass is fixed to the fiducial value 
fu = 0. 
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Figure 11. The marginalized error ellipses in a sub-space of (u>o, fv), where both the parameters are included in the Fisher analysis. Solid ellipse represents 
the error when the FoG shape is completely known (i.e. "FoG shape known"), while the other ellipses show the errors when the FoG effect is approximated 
with Gaussian or Lorentzian forms "with offset" (dotted and dashed respectively). The imperfect modeling of the FoG effect biases the results different from 
the input values denoted by cross symbols. The lensing information is added for all plots. The different panels differ in the hypothetical spectroscopic survey 
(BOSS or PFS) and the maximum wavenumber fc ma x used in the Fisher analysis (C = 0.2 or 0.7). 
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Figure 12. Similar to the previous figure, but for the joint experiments of HSC+BOSS+PFS+Planck (upper panels) and Euclid+Planck (lower). 



5 SUMMARY AND DISCUSSION 

Our lack of understanding of the relationship between the positions of galaxies in redshift space and the underlying distribution of matter is 
the largest systematic source of errors in the analysis of redshift surveys. If the dominant galaxy (e.g., the brightest cluster galaxy) in each 
halo was at rest in the center of the halo, then there would be minimal Finger-of-God (FoG) effects in a properly selected sample, nor any 
suppression of the power spectrum amplitude due to finite size of galaxy halos. We have quantified the deviations from this simple universe 
in terms of galaxy radial distribution, the probability of finding a galaxy at a given distance from the center of the halo. Since galaxies not in 
the centers of the gravitational potential are moving relative to the potential, thus the amplitude of this displacement is directly related to the 
amplitude of the FoG effect. 

In this paper, we have showed how galaxy-galaxy lensing observations, which measure the relationship between galaxy positions and 
the distribution of dark matter in a statistical sense, can be used to measure this profile and reduce the uncertainties in analyses of galaxy 
power spectrum measurements. We have studied a case of luminous red galaxies (LRG), especially the dominant LRG in each halo (DLRG), 
as a working example for demonstrati ng an expected p erformance of our proposing metho d. We have assumed DLRG off-centered profiles 
motivated by the weak lensing study ( Og uri et aljEoioh and the analysis (Ho et al. 2009) studying offsets of LRGs relative to the peak in 



X-ray brightness. Our DLRG radial profile model predicts that the FoG effect suppresses the DLRG power spectrum amplitudes at k ~ 0.1 
and 0.2h,/Mpc by 10 and 30%, respectively, much greater than the statistical precision of power spectrum measurements for ongoing and 
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upcoming spectroscopic surveys (see Fig. 8). The halos with masses ~ 10 Mq make the dominant contribution to the FoG effect on the 
DLRG power spectrum. 

If we can measure the radial distribution, then we can model this effect and remove (or at least significantly reduce) this source of 
systematic uncertainty. We have illustrated this correction for ongoing and upcoming surveys: the BOSS survey, the Subaru imaging (HSC) 
and spectroscopic (PFS) surveys and ESA's Euclid mission (see Table [TJl. For example, the Subaru HSC imaging survey covering an area 
of 2000 square degrees can accurately measure the DLRG-galaxy lensing profile down to scales of an arcminute enabling an accurate 
characterization of the DLRG off-centered profile (see Fig.|6]and Table|2j- As we have showed, this radial distribution can be used to correct 
the FoG effect and significantly improve parameter estimation from the redshift-space power spectrum measurements. We have illustrated 
this effect by discussing the impact of FoG effect on the parameters wq and /„ as well as the growth rate f B = d In D/d In a at each redshift 
slice. We have found that combining the lensing measurements with the redshift-space power spectrum can improve the constraints on f g up 
to a factor of 2, compared to the case that the FoG effect is considered unknown, if the power spectrum information down to fc max ~ 0.1 
or 0.2/i/Mpc is included (see Fig[9jl. Including the lensing correction also improves measurements of wo and /„ by 5-25%. Such a factor of 
2 improvement in the parameter corresponds to a factor of 4 larger survey volume, and therefore the combined imaging and spectroscopic 
surveys can bring a huge beneficial synergy. 

The analyses in this paper have made a number of simplifying assumptions that would have to be improved in a realistic analysis of a 
large spectroscopic survey. We have assumed a linear model for redshift distortion effect on halo power spectrum (the Kaiser formulation 
in Eq. (5)). The qual ity of the rapidly improving redshift data requires more refined model bas ed on extended perturbation theory as well as 
simulations i Scoccimarro 2004 ; Taruva et alJboiol : iTang et aT 2011 ; Sato & Matsubara] 2011 ). We have also assumed a linear bias model. 
Saito et all J201 1 ) found that non-linear bias corrections are important for the SDSS DLRGs. 



In this paper, for simplicity, we have assumed that the velocity dispersion of DLRGs is determined by an enclosed mass (mostly dark 
matter) at a given radius and that the velocity distribution of DLRGs obeys a Gaussian distribution in a statistical average sense (see Eqs. 1141 
and 1171 ). However, the velocity distribution should be by nature sensitive to various complicated physics and galaxy formation physics, 
so it is still difficult to understand based on theoretical studies. Therefore an observational approach to tackling this issue is rather more 
adequate. For example, one may be able to use spectroscopic redshifts of galaxies, from a survey data itself or from a dedicate d survey with 



optimized fiber position i ng, to explore the v elocity structures of DLRGs around massive halos based on the stacking analysis dBecker et al 



20071: IwhiteetalJboiol : Skibb a et alj|201 11) (also see Lam et al. in preparation). Or an alternative approach on the analysis side is to use 



an empirical model of the FoG effect which is given by a multiplicative function of the Gaussian FoG form and a perturbative functional 
form [1 + ao(kfi) 2 + ai(kfi) 4 +•••]> where the coefficients ao, ai are treated as nuisance parameters. The model defined in this way can 
have more degrees of freedom to describe possible complicated scale-dependences in genuine FoG effects 1 Tang et all 2011 , for a similar 
approach). Then, by using the generalized FoG form, we may be able to derive more robust, unbiased cosmological constraints including 
marginalization over the nuisance parameters, but still using the lensing information of off-centered DLRGs to constrain the Gaussian FoG 
part as we did in this paper. 

There is more work needed to explore the effects of DLRG offsets on galaxy power spectra and lensing. Since galaxies of a given type 
lies in a range of halo masses, it is important to trace the dependence of the DLRG clustering and the FoG effect on halo mass scale: for the 
clustering, halos with typical mass scale of ~ 1O 13 M gives a dominant contribution, while the FoG effect arises mainly from more massive 
halos with masses around 1O 14 M0. It would be useful to select a sample of massive halos focused on the high mass range to more accurately 
explore the offset of the DLRG from the halo centers as defined by optical, X-ray and SZ data. In this approach, weak lensing can play an 
important role, as various halo cent er indicators can be used to monitor weak lensing signals at small angular scales in both individual cluste r 
lensing and stacked analysis bases 1 Johnston et alj2007 ; Leauthaud et al .l20ld l Okab e et al ,l2010l :l lOguri et al ■l2010l:IOguri & Takadalboill) . 

Besides the cosmological use, the off-centering information of DLRGs should be very useful to develop a more physical understanding 
of the dynamical processes of DLRGs within the main halo. DLRGs are likely to reside on one of the most massive sub-halos within the 
main host halo with Si 1O 13 M0. Therefore DLRG within each halo tends to sink towards the halo center due to dynamical friction, which is 
one of basic explanations for a more centrally concentrated distribution of DLRGs within halos compared to the dark matter distribution. The 
dynamical processes and assembly histories of DLRGs, within a given time scale of cluster-scale halos, are a key information to understanding 
the nature of DLRGs in the context of CDM dominated structure formation scenario. Again the galaxy-galaxy lensing can offer a new method 
of tackling these issues. 
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APPENDIX A: ISOTHERMAL VELOCITY DISPERSION MODEL 



In this appendix, we show an alternative model of the velocity dispersion of DLRGs following the theory in iBinnev & Trema ine (2008), 
which differs from the model in Section l2.3.3l 

Let us begin our discussion with assuming that the phase s pace density of DLRGs, which reside in host halo of mass M, obeys an 
isothermal distribution (Eq. [4-116] in lBinnev & Tremaine! ( 20081) ): 



/(r,v;M) 



Pi 



(2^ LRQ (r;M))3/2 



exp 



*(r) 



72 



,(r;M) 



(Al) 
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where $(r) is the gravitational potential determined by dark matter distribution, c?\ i ao (r; M) is the ID velocity dispersion of DLRGs and pi 
is the normalization constant. Here we assume a spherically symmetric distribution for the DLRG distribution in a statistical average sense. 
Although DLRGs obey such an isothermal distribution in the potential well of dark matter, this is an alternative model to estimate the DLRG 
velocity dispersion, so let us continue our discussion. 

Integrating the equation above over velocities yields the radial profile of DLRGs, which should be equivalent to the off-centered profile 
in our language: 

( *(V) \ 

p off (r; M) = pi exp -= )-!— . (A2) 



This equation gives us an inverse problem: once the off-centered profile and the gravitational potential are given, the velocity dispersion 
0ti,isois determined. 

Assuming an NFW profile for the total mass profile within the halo, the gravitational potential is given via Poisson's equation as 
Id/ 2 d*(r) 



r»dr V r ^)=-^G PD u(r), (A3) 
Hence the potential is given as 

, GM(< r') 



*(r) = - dr 



2* 2 v (r;M) /(r/ ;; ) / , (A4) 



m nfw (c)/c 

with 

m n f w (a;) = ln(l + x) 



1 + x' 



1{x) s / dx 'Hh^l = i _ Ml+x), (A5) 



Here we have rewritten the equation above in terms of the virial velocity dispersion a v (r; M) defined in Eq. il4l . 
Inserting Eq. &A4\ into Eq. dA2t yields 

t *n ( 2ag(r;M) I(r/r 3 ) \ 
p ff(r; M) = pi exp ^ --— . (A6) 



J v,i so ( r ; M ) m nfw (c)/c 

The normalization constant pi is determined by imposing that the velocity dispersion a V: i so is the same as a v at the limit r — > 0: 

^ Poff(r pi v)exp(2) (r piv /r s < 1) (A7) 

where r p i v is a pivot radius we choose to impose the conditions that r p i v is very small as well as a„,i so = a v at r — r p i v . Hence, the velocity 
dispersion for an isothermal distribution model can be given in terms of the off-centered profile as 

^ iso (r; M) = ^(M)-^L- ""f'^ ) } • CA8) 

m nfw (c)/c [ 2 \Pos{r P i v ; M) J _ 

Fig. lAll shows the velocity dispersion for this isothermal distribution model, comparing with our fiducial model based on a simple virial 
theorem (Eq. 11141 ). To compute this velocity profile, /r s = 10 4 . However note that the dependence of r p i v is small because 

the dependence is logarithmic as implied in Eq. dA8t . The figure shows that the isothermal model gives a smaller velocity dispersion than 
the fiducial model, and therefore the resulting FoG effect is smaller than the results shown in the main text, by about factor of 2 in the power 
spectrum amplitudes. 
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Figure Al. Comparing the models of DLRG velocity dispersion within a halo with 10 14 Mq and at redshift z = 0.45. The two models to be compared are a 
model based on the virial theorem (Eq. 1141 ) and a model derived assuming an isothermal velocity distribution given the off-centered profile of DLRGs and 
the dark matter mass profile. The left- and right panels show the results assuming the Gaussian and NFW off-centered profiles, respectively. The model given 
in Appendix lAl gives a smaller velocity dispersion than our fiducial model. 
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